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Abstract 

In this paper we study the effect of next-nearest-neighbor hopping on the 
dynamics of a single hole in an antiferromagnetic (Neel) background. In the 
framework of large dimensions the Green function of a hole can be obtained 
exactly. The exact density of states of a hole is thus calculated in large di- 
mensions and on a Bethe lattice with large coordination number. We suggest 
a physically motivated generalization to finite dimensions (e.g., 2 and 3). In 
d = 2 we present also the momentum dependent spectral function. With 
varying degree, depending on the underlying lattice involved, the discrete 
spectrum for holes is replaced by a continuum background and a few reso- 
nances at the low energy end. The latter are the remanents of the bound 
states of the t — J model. Their behavior is still largely governed by the pa- 
rameters t and J. The continuum excitations are more sensitive to the energy 
scales t and ti. 
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I. INTRODUCTION 



The dynamics of a hole in an antiferromagnetic background is a physically interesting 
but complex problem with a long history in particular in the theory of magnetic 

semiconductors 0. In solid ^He it precisely corresponds to the motion of vacancies at 
low temperatures. Due to the fact that in lightly doped High-Tc superconducting materials 
holes move in the antiferromagnetic background of the original undoped insulator the 
problem received wide attention again, leading to exciting new physical insight. Nevertheless 
a complete solution still does not exist. The theory has recently been reviewed by Yu et al. 

i- 

The physical problem is most easily understood for an Ising antiferromagnetic (i.e. Neel) 
background with exchange coupling J in two dimensions. As the hole moves, it breaks spin 
bonds. A straight-line motion over a distance x thus leads to a string of broken bonds 
and thus to an increase in energy proportional to Jx. At the same time strings unwind 
themselves if the hole moves around a loop one and a half times @]. Nevertheless, since the 
latter is a rather high order process in the hopping amplitude t, a hole has been considered 
to be bound in a linear potential 0J^. The resulting hole spectrum consists of a sequence 
of discrete eigenvalues. For a t — J model these can be seen to scale with t[J/tf'^^ from 
essentially dimensional arguments. Inclusion of Heisenberg-like spin flip interactions leads 
in turn to delocalization of the hole and the formation of an effective quasiparticle band 
whose width is of the order of J. This scenario has been studied extensively in one and two 
dimensions using various analytic and numerical techniques |1T0|-|T6i[7| . 

We shall argue below that physical considerations suggest an approach that involves an 
additional energy scale ti, the hopping amplitude to next-nearest-neighbor sites. A realistic 
band structure simply cannot be fit by an energy dispersion that originates solely from 
nearest-neighbor hopping |T7|-[19|. In fact, estimates for \ti/t\ range from approximately 
0.15 for La2_xSr^Cu04 to roughly 0.45 for YBa2Cu307_2^^ [|17|, indicating that ti is far from 
being a small energy scale. The significance of a ti term is profound. In a bipartite lattice, 
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motion via ti leaves the spin configuration invariant. Hence the particle becomes free to move 
without having to distort the spin background. While Trugman loops ^ can be viewed as 
effective ti processes where ti ~ t^, we expect that a bare ti exists as an intermediate energy 



scale. We are thus led to investigate the t — ti — J model and to ask the question pO[]:/tow 
does ti alter the physical picture described ahovel 

For the t — J model the problem of a single hole in an antiferromagnetic background 



was recently solved in the limit of large dimensions, — > oo by Metzner et al. ||22 



for J = and by Strack and VoUhardt [^,^] for J > 0. They were able to show that 



string unwindings, whether because of Heisenberg-like spin flip interactions or as a result of 
Trugman loop motion, become negligible in the limit of infinite d. Hence approximations 
such as the Brinkman-Rice retraceable path approximation (rpa) or the linear potential 
approach [^,0 become exact at c? = oo. As a result, the density of states for a hole and 
its dynamical conductivity can be calculated exactly for d = oo. Moreover, systematic 
improvements may be obtained through an expansion in powers of l/rf p2|-p^. They show, 
for example, that the results of the rpa in finite dimensions correspond to a self-consistent 
l/(i-expansion of the exact d = oo result pB| , |2^ . Furthermore, starting from the rpa and 
resumming the Nagaoka expansion in terms of non-retraceable skeleton paths dressed by 
retraceable-path insertions Miiller-Hartmann and Ventura recently obtained an almost 



quantitative solution of the one-particle problem in all dimensions. There is good agreement 
with the numerical results by Zhong et al. [^, which were obtained for a double chain and 



a square lattice with a new Lanczos iteration scheme. 

For the t — J model in infinite dimensions the Green function is local and satisfies a 



self-consistent equation [23 



which may be solved explicitly in terms of Bessel-functions Here \i) represents the Neel 
state with a hole at site i. The parameters t* and J* are appropriately scaled energies whose 
precise definition will be given below. Eq. (|l]) coincides with the approximation by Kane et 



al. [T^l; the latter is therefore exact in large dimensions. 

In the light of its success in capturing several of the essential physical features of hole 
motion in the t — J model, we employ the limit of large dimensions also to the t — ti — J 
model. Exact expressions are derived for the hole Green function, Gij{z), at d = oo and 
a natural extension to finite dimensions is provided. The latter continues to possess all 
required analytical properties, including spectral sum rules. 

The paper is organized as follows: Sec. II contains a description of the Hamiltonian 
and the necessary parameter scalings for d = oo. In sec. Ill, we derive the Green function 
for holes on a hypercubic lattice in ci = cxd, for a Bethe lattice with Z = oo, and for a 
generalization of our results to finite d (= 2, 3). Sec. IV contains a discussion of our results, 
followed by a summary in sec. V. In an appendix we address the analytical properties of 
our finite d generalization. 



II. THE T-Ti-J MODEL AND THE LIMIT OF LARGE DIMENSIONS 

The Hamiltonian we consider, written for convenience in terms of creation and annihila- 
tion operators for holes, is given by 

= nt + nt,+nj. (2) 

Here /ijo- (^L) annihilates (creates) a hole with spin a on site i , and Si is the spin at 
site i. The operators n^^ = hl^hi^ , which appear in the hopping terms, ensure that no 
electronic double occupancy occurs in the course of hopping. As always, denotes 
nearest-neighboring sites, while the notation represents next-nearest-neighbors. 

We adopt Eq. (§) as our basic model, whose microscopic justification traces back to the 
Hubbard model with a next-nearest-neighbor hopping term added. Upon carrying out the 
standard canonical transformation for a large U and low hole dopings ||28|, a number of 
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terms are generated. In addition to the J-term (J = At^ /U) specified above, we recover 
an antiferromagnetic coupling Ji = Atl/U between next-nearest-neighbor spins, as well as 
a number of three-site terms. The latter include also processes where t and ti hoppings 
are combined. Clearly, a strong enough Ji coupling will lead to frustration in the low 
temperature antiferromagnetic order at half filling. Yet, by focusing on the physical range 
where t constitutes an energy scale larger than ti , the role of Ji is mainly to renormalize 
the effective coupling constant J |2^. We therefore incorporate these effects in Eq. (|^) into 



the (independent) parameter J. At the same time, for large U, ti itself may very well be 
appreciably larger than J. 

Figure 1 depicts the hole trajectories in the antiferromagnetic background which become 
important in the presence of a ti term. Here we are considering paths that bring the hole 
back to its starting position without altering the spin configuration. The most significant 
contribution arises from paths shown schematically in Fig. 1(a), since they leave the spin 
configuration invariant. At the same time a new type of loop motion, sketched in Fig. 1(b), 
appears. Here, by combining two t and one ti hoppings, a triangular loop motion is sufficient 
to restore the antiferromagnetic background. This must be compared with the case for 
ti = 0, where three circulations around a plaquette are necessary for the ordinary Trugman 
loops. 

Next we introduce the limit of large dimensions. At large d the model parameters are 



scaled so that [21-23 



t = ^= , ti = 4 , J = — . (3) 
V2d d ' 2d ^ ^ 



Keeping t*,tl and J* fixed, a meaningful d ^ oo limit emerges. One consequence of such a 
limit is that spin flip processes become negligible, so that spins may be regarded as being 
Ising-like [^. In addition, without ti, only retraceable paths are allowed up to 0(1/ d'^) 
P2|| . Now, however, both paths displayed in Fig. 1 yield non- vanishing contributions since, 
for example, hopping on a triangular loop (see Fig. 1(b) ) is proportional to t^ti ~ l/c?^, 
while the number of different embeddings for this path is proportional to d"^, such that the 



product of the two remains of order unity even when d approaches infinity. More generally, 
by an appropriate combination of t- and ti-processes (or using only ti-processes) the hole 
can now move on closed loops of arbitrary length even in = oo. Obviously the resulting 
motion of the hole, consisting of simple, retraceable paths and ti-generated, background 
restoring loops, can be very complex. 



III. THE GREEN FUNCTION FOR HOLES 

The hole Green function is given by 

G^,{z) = (Aj^lj) . (4) 

Clearly i and j in Eq. must lie on the same spin sub-lattice (taken hereafter to be the 
spin up sub-lattice), otherwise Gij trivially vanishes. 

Typically, Gij{z) is calculated using Nagaoka's path formalism In this formulation 
[z — 7i)~^ is expanded in powers of Ti/z, and the different contributions are identified with 
all the different background restoring paths that extend from site j to site i. Later on we will 
indeed return to this approach, but first we wish to account for the fact that motion via ti 
alone is completely unrestricted. For this purpose it is useful to construct two complementary 
projection operators: 

^ = EI^)(^I > Q = l-V. (5) 

Here the summation over i is confined to the spin up sub-lattice. The hole Green function, 
(I), corresponds then to the operator 



The latter may be expressed, with the aid of a familiar operator identity ||30 



as 



G{z) = ^r- . , (7) 

G(o)-Hz)-nz) 



where 



G^^)[z)=v--—v , ±{z) = n„-——n,, (8) 

and 

Hpp = vnv , ngg = QnQ , etc. (9) 

Here G^^\z) accounts for all pure ti processes, while all processes involving t are incorporated 
in the form of a self-energy. This separation into spin ground-state processes and those 
involving spin excitations is very useful in simplifying the algebra. For example, Hpg and Hqp 
in the numerator of £(2;) necessarily require only Tit and Tij^. However in the denominator, 
Hqg contains all processes where the hole moves in an excited spin state, including motion 
via ti. 



A. The self-energy at large dimensions 

Physically, all the complexity of the hole dynamics is contained within the self-energy. 
Yet at finite dimensions (i.e., d = 2,3), a complete solution for Sjj(z) remains out of reach. 
We show, by contrast, that the problem can be solved exactly at infinite dimensions, due to 
a sequence of simplifications. Primarily the direct effect of spin flips on the motion of the 
hole vanishes as c? — 00 p2| , p3| and hence the self-energies can be rewritten as 



{i,k) z- c;[nt + nt, + njjc; 



Next, as usual in d = 00, Tiij{z) becomes site-diagonal |2T|. This characteristic large 



dimensions feature follows from the application of Nagaoka's path formalism to (plQl). A typ- 
ical path describing ^ij{z) begins with an excitation of the ground-state spin configuration 
at site j, and its restoration only at the final step when the hole reaches site i. In between, 
the spin background is continuously excited during the entire hole motion. Consider now the 
case where i and j are two different sites (Fig. 2). Then the fact that the spin background 
is eventually restored at a site i different than j demands that the latter must be revisited 
by the traveling hole at some intermediate step. Otherwise, the down spin it acquired in 
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the initial step will never be mended, as it should, to an up spin. Since the overall spin 
background remains excited through all intermediate stages, this implies that some other 
site / on the up sub-lattice is forced to acquire a down spin, in between the two hole visits 
to site j. Moreover, after the spin at site j is finally corrected to an up one, site / has 
to be revisited in the course of restoring the ground-state spin configuration. A schematic 
description of the path is shown in Fig. 2. It consists, among other things, of at least three 
segments where the hole propagates between the initial site j and the intermediate site /. 
From the standpoint of large dimensions this is two segments too many. Thus the effect 
of the spin background on the off-diagonal components of the self-energy is to reduce their 
contribution by an integer power of l/d"^. Indeed the path illustrated in Fig. 2 does not have 
the topology of a "loop tree" ||2^ and hence does not contribute at d = oo. 



The third and final simplification at large dimensions occurs in the calculation of the 
remaining on-site self-energy. In order to evaluate Tjii{z) we focus on the following functions 

= -e(«.+k + H,.)e ">^'i'> 

{k, I being nearest neighbors of i), and compare them with 

Gkiiz) = (Neel|4^ J /^u|Neel) . (12) 

z — rit -r rLti 

Here |Neel) designates the antiferromagnetic spin ground-state with no holes present. Eq. 
(|l^) is no other than the hole Green function, only written in terms of a hole in the down 
spin sub-lattice rather than the up spin one. Ti\^'^\z) and Gki{z) differ in two respects: (1) 
the Q projection operator appears in the denominator of Eq. ([TT|), and (2) the spin at site i 



is flipped in Eq. (|TT]) to a down spin. Both differences are conveniently accounted for in the 
limit of large dimensions by resorting once again to Nagaoka's path formalism, and pursuing 
the same line of reasoning as before. Here one can see that all possible hole trajectories 
where site i is involved are reduced by an integer power of l/d for both and Gki{z). 

Hence the spin at site i may be regarded as being frozen in each of the two cases. Such an 
observation yields profound consequences. To begin with, the down spin at site i in the case 

9 



of s|f''-*(z) ensures that the spin background remains constantly excited. Therefore the Q 
projection operators entering Eq. ( pA]) are automatically satisfied. In addition, since apart 
from the flipped spin at site i the initial and final spin configurations in Eq. (pA]) are identical 
to those of Eq. ([T2|), we precisely recover the same path contributions for both cases. Thus 
the only distinction between the two cases exists in their Ising energy, which happens to be 
J*/2 larger for T.[i'^\z). This is a result of the broken bonds between the flipped spin at site 
i and its surrounding nearest neighbor spins. All in all we arrive at the following relation, 
which becomes exact for large dimensions: 

Er)(^)=G.K^-y). (13) 
Finally, the combination of Eqs. ([T0|) and (|l^) yields 

m^^d^)=t'j:j:GM{z-^). (w) 

(i,k) {i,l) 

Note that no restrictions on the nature of the underlying lattice were involved in the 
derivation of equation ([1^), hence it equally applies to all bipartite lattices in the limit of 
infinite number of nearest neighbors. 

B. The hole Green function at large dimensions - hypercubic lattice 

After establishing that the self-energy is site-diagonal, Eq. is readily solved for a 
hypercubic lattice by Fourier transformation: 

with 

= cos(A;„) cos{km) ■ (16) 

Here k lies within a reduced Brillouin zone, and the lattice constant is taken to be unity. 
The spatial Green functions are obtained from the inverse transform of Eq. (plsl). It yields 
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G,,{z) = Gif{z-E{z)), (17) 

where G^-^^ is the free, tight-binding Green function due to next-nearest-neighbor hopping. 
Next we make use of the cubic symmetry, and express ^{z) from Eq. ([T^) as 

E{z) = [g{z - y) + 2{d - l)G,{z -^-^)+ G,{z (18) 

Here the notation G = Ga and Gi was introduced for the on-site and next-nearest-neighbor 
Green functions, respectively, while G2 corresponds to the case where the two lattice sites 
lie on the same axis, yet with a separation of two lattice constants (i.e., third-order- nearest 
neighbors). The former two functions are further related through an identity that links the 
corresponding local and next nearest neighbor tight binding Green functions. Indeed for 
any dimension 

2(.-l)GrW^£q£l^. (19) 

where Gf = G^^\ Thus, since G2{z) vanishes at large dimensions as 1/d, Eqs. ([T8| ) and 
(O) reduce at the limit of c/ 00 to 



Together, equations (|T^) and (|20|) combine to give an iterative equation for the local 
Green function G{z), which is exact a.t d = 00. It must be supplemented, however, with the 
corresponding d = 00 expression for G^^^ (z) which may be obtained (using the corresponding 
density of states E^) as P2I 



G<».(.) = Zl , ^\ ,^ W J^] ^ (21) 



2«;^(z + (;)/2t; VV 2«; 

Here W{z) is the scaled complementary error function of a complex argument |3^, and the 
sign of the square root is chosen so that it lies in the upper half plane. 

C. The hole Green function on a Bethe lattice 

In case of a Bethe lattice, the definitions of the scaled parameters are modified as 
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where Z is the coordination number. Once these adjustments are carried out, equations (|T 
and ( pUD continue to hold as the exact Z = oo equations for the hole Green function Gij{z) 
where now, however, G^'^\z) in Eq. (^) is given by [§] 



= ■ (23) 



D. Extension to finite dimensions 

A systematic 1/d expansion about the d = oo limit constitutes in general a formidable 
task, and in that respect the problem of hole motion in the t — ti — J model is no exception. 
Here we consider a hypercubic lattice and suggest an alternative, approximate way to bridge 
between d = oo and finite d. 

We base our approach on the same equations that were derived for large dimensions, yet 
by two adjustments. First, the G[f{z) Green functions are replaced by their corresponding 
(i-dimensional expressions. Second, we retain the contribution of G2 to the local self-energy, 
i.e. use Eq. (|18|) rather than Eq. (|20|) . Though the latter contribution completely vanishes as 
(i — s> 00, it is essential to keep it at finite d in order to avoid a non-physical spectral function 
(note that for < 00 a shift occurs in the band edge of the spectral part of G^^\uj — i5)] see, 
e.g.. Fig. 3). We elaborate on this point in the appendix. The complete set of equations 
employed for a finite d consists, therefore, of Eq. (|l^, together with the following version 
of Eq. (m-. 



S(.) = { Giz - + + (24) 



We further note that the only input function required at d = cxd is G^^\z), while for a 
finite d it is joined by G2'\z). For d = 2,3, both functions can be computed using complete 
elliptic integrals of the first and second kind, provided they are properly continued to the 
entire complex plane |]3^. In the appendix we show that all analytical properties required 



from Gij{z) are indeed conserved within the approximation scheme suggested above. 
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IV. NUMERICAL RESULTS 



A. The bare density of states 

Figure 3 shows the bare density of states (DOS) due to pure ti-hopping on a hypercubic 
lattice in d = 2, corresponding to the DOS for a square lattice, in = 3 and d = oo [see Eq. 
(21)], corresponding to the DOS for a fcc-lattice in those dimensions, and, finally, for a Bethe 
lattice with Z = oo, Eq. (23). There is a divergence at the lower band edge (for tl > 0), 
either with a square root singularity for d = oo and the Bethe lattice, or in a logarithmic 



manner 



3^ for (i = 3. For d = 2 the logarithmic singularity is at the band center. The 
curves for d = oo and Bethe lattice are clearly very similar, although the bandwidth for a 
Bethe lattice remains finite. In contrast, d = oo possesses an exponential tail that extends 



to higher energies For d = 2 and 3 there is a slight shift in the lower band edge as 

mentioned earlier in sec. III.D. Except for d = 2 the cases look qualitatively similar. For 
tl < the bare density of states are inverted (i.e., u —u). The divergences then move to 
the upper band edge, with the exception of the d = 2 case which is symmetric in u. 



B. Infinite dimensions 



Fig. 4 shows the effect of an increase in > on the DOS for the motion of a hole in 
the infinite-dimensional t — ti — J model at fixed J* /t* = 0.4. For tl/t* = 0.05 we see only 
little changes from the case where = 0: there are new satellites in-between the original 
bound state energies, starting a.t ljq + J* /2, where Uq is the lowest bound state energy; the 
sharp delta peaks for tl = have now acquired small width. (Note, however, that we expect 
the realistic energy hierarchy to be t* > tl^ J*.) At tl/t* = 0.2 a continuous spectrum has 
developed. The sharp structures near the lower band edge are derived from the divergence 
of the bare DOS. At tl/t* = 0.4 only the lowest resonances appear discernible. The rest has 
been absorbed into the incoherent continuum of excitations. Any further increase in tl leads 
to additional decrease of the resonant peaks and a further development of the underlying 
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continuum. It looks as if one is approaching the general shape of the bare (continuous) 
density of states. 

For tl < 0, Fig. 5, there is a dramatic change in the DOS. By tl/t* = —0.2 only weak 
traces are left of the bound states at tl = 0. The DOS is more symmetric and less structured 
compared to > 0. It corresponds to an incoherent band with an exponentially small tail 
that extends to negative energies. The origin of the strong effect that a negative has on 
the hole motion lies in the unboundedness of the bare DOS for uj — oo. As a result, the 
self-energy acquires a nonzero imaginary part throughout the negative energy region. 

Figure 6 displays the DOS for a Bethe lattice with Z = oo for J*/t* = 0.4 and tl/t* = 
±0.4. Again, well-developed incoherent continua are recovered, similar to the ones obtained 
for d = oo. The resonances persist, however, even for tl/t* = —0.4. They seem to be 
the dominant structures for smaller values of {\tl\ ^ 0.2). Notice also the shift in the 
threshold energy between tl > and tl < 0. For tl < the resonances shift to lower energies 
with increasing while in case of > the lower edge is quite insensitive to variations 
in tl. The difference between d = oo and the Bethe lattice is mostly due to the finite width 
of the bare spectrum of the latter. 

C. Finite dimensions 

Next we turn to finite dimensions. Figure 7 depicts our generalization to d = 3 for 
J*/t* = 0.4 and tl/t* = ±0.2, ±0.4. The width and magnitude of the incoherent spectra 
seem to be similar for t* < and tl > 0, as for all lattices discussed here. The lower edge 
is found, once again, to have no significant shift with t* for tl > 0, while it decreases with 
decreasing tl for tl < 0. The qualitative behavior of the DOS for tl/t* = 0.4 is similar to 
that obtained by Bala et al. for d = 2 using a very different method. The effect of a 
smaller J*, corresponding to what might be expected for realistic material parameters, is 
shown in Fig. 8; for J*/t* = 0.1 and tl/t* = —0.4, we recover a broad incoherent spectrum. 
At the same time, resonances can be clearly seen at the lower edge, with widths that are 
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smaller than the separation between adjacent peaks. The low-energy features thus persist 
even for J* appreciably smaller than t\. 

We now examine more closely the low-energy spectrum for two dimensions. For d = 2, 
the general features described above are seen as well. The effect of increasing \tl\ is to 
broaden the original bound states into resonances, eventually leading to continuum in the 
DOS. There are, however, two major differences. The first is that a larger is needed 
to produce qualitatively similar results in c? = 2, as compared to d = 3. The reason is 
simple: the overall bandwidth of the bare DOS is smaller for = 2 (see Fig. 3). The other 
difference is that in d = 3 a negative tl is seen to be more effective in smearing out the 
residual details of the bound states. In d = 2, the roles are slightly inverted. This has to do 
with the symmetric character of the bare DOS in c? = 2 as opposed to the strong asymmetry 
for d>3. 

We have also studied the fc-dependent spectral function A{k, u) for d = 2. Figure 9 
shows A{k,uj) for k = (|, |) and J*/t* = 0.2, tl/t* = -0.6. It looks quite similar to the 
spectral function calculated numerically by Liu and Manousakis [0 for the t — J model on 
a finite size lattice {J* /t* = 0.2 corresponds, for = 2, to J/t = 0.1). The lowest peak 
in this curve corresponds to a coherent state. Indeed all fc-states within a margin of J*/2 
from the lower band edge have vanishing imaginary part of the self-energy. Physically, this 
corresponds to their inability to decay by emitting or absorbing spin excitations, which is a 
consequence of the Ising nature of the spins. For tl < the dispersion of this coherent band 
is similar to that of the bare energy in Eq. (p!6D, with a strong renormalization due to J*. 
To be more specific, we find that 

-> J* -> 

E{k) = —f (k; tl) + constant . (25) 



2 

For a fixed tl and different J*, the eigenvalues scale with J*/2 in that E{k) — E{ko) / (^^) 
all lie on the same curve. This is demonstrated in Fig. 10(a). The function /, however, 
differs from Eq. (|TB|) in that the upper band edge is deformed towards lower energies for 
nonnegligible values of t* [see Fig. 10(b)]. For small values of tl, though, / is essentially 
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linear in t*. The overall integrated spectral weight of the entire band can be expressed as: 
where p^^^ is the bare DOS for d = 2. Following E{k), Eq. (^) also scales with J*. One con- 

— * 

sequence of the self-energy being /c-independent in our approximation is that the minimum 
of E{k) always occurs at the F point, k = (0, 0). 

For positive t* in c? = 2, one notable feature is that for moderate t* {tl ^ J*), when the 
low energy resonance in the DOS has almost become part of the continuum, there are wave 
vectors for which coherent state solutions no longer exist. Three types of solutions are then 
recovered. At low frequencies, and over a frequency range of Auj = J*/2, true coherent state 
solutions are found. These are centered around the band minima at k = (vr, 0). At slightly 
higher frequencies, the imaginary part of the self-energy becomes nonzero in a quadratic 
manner. The solutions are thus quasiparticle-like. In the neighborhood of A; = (0, 0) — 
the upper band edge — and over an area in fc-space that increases with increasing t^, 
the solutions are removed from the low energy side and are shifted to the incoherent high 
frequency side of the spectrum. 



V. DISCUSSION 

In this paper we explored the consequences of a next-nearest-neighbor hopping ti on 
the motion of a single hole in a quantum antiferromagnet. The system is described by a 
t — ti — J model. Our motivation stems both from the need of a ti term to characterize 
the observed electronic structure |ll7Hl9[] , and from the growing numerical evidence that 



the effective quasiparticle dispersion generated by the t — J model can be parameterized 
by means of effective next-nearest-neighbor and third-order-nearest-neighbor hops. Our 
approach is based on the limit of large dimensions, where the problem is solved exactly. For 
d = oo, ti survives as the only mechanism that delocalizes the hole. The principal effect of 
the introduction of ti is thus seen in the following steps: (a) For tl = (and J* 7^ 0), we have 
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only bound states that are typically separated by t{J/ty/^; (b) as t* increases (t* ^ J*), all 
peaks in the density of states acquire width, and satellites appear at intermediate energies 

n,mJ* /'^ (here Un is the energy of the ri-th bound state and (in,m ~ is^ 
t\ further increases {t\ ~ J*), the satellites and resonant peaks merge, and a continuum is 
formed. Hence ti restores the incoherent Brinkman-Rice |^ spectra at high energies, while 
generating a quasiparticle dispersion at low energies. This is most clearly seen within our 
analysis in ci = 2. For \t\\ ^ J* and uj close to the low-energy edge, we may replace the 
self-energy in Eq. ( |T^ ) with its t\ = Q limit. Then, using a dominant pole approximation 

m ' 

+ = ^, . (27) 

where ujq and are the = ground-state energy and quasiparticle amplitude, respectively. 
Hence a quasiparticle band immediately appears, with a characteristic renormalization Zq. 
Its bandwidth is given by ^ZqVI (the factor 4 is for d = 2), while the /c-independent effective 
mass is proportional to 1/ (Zot*). The effect of J* enters via the quasiparticle amplitude 
Zq. In the intermediate and large regimes for t\/J*, the bandwidth is practically governed 
by J* (see, e.g.. Fig. 10). For t\ < 0, the entire Brillouin Zone is included in this band, 
and the quasiparticle dispersion scales with J* . At moderate values of t\ > 0, quasiparticles 
solutions are removed from the vicinity of the F point. In both cases, the frequency range of 
width J*/2 at the lowest band edge accommodates only coherent states. This corresponds 
to the inability of holes to loose energy by flipping the spins. The latter is a consequence of 
large dimensions where the original Heisenberg spins become Ising-like, i.e. quantum spin 
fluctuations are suppressed. Quantum fluctuations are known to have two important effects 
in the t — J model. Together with Trugman loops they act to delocalize the hole, giving 
rise to an effective quasiparticle energy dispersion that can be fitted by effective hops to 
next-nearest-neighbor and third-order-nearest-neighbor sites. Quantum fluctuations are also 
responsible for the quasiparticle-band minimum at k = (f , |) (see, e.g., Zhou and Schultz 
||T6| for a clear demonstration of this point). The introduction of a ti-term to the model, 
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motivated by physical considerations, simulates quite well the dispersion, i.e. delocalization, 
of the hole but cannot produce a band minimum at k = (f , f). 

In summary, we have provided a simple analytic theory for the hole-motion in an anti- 
ferromagnetic background in the presence of next-nearest-neighbor hopping. The theory is 
exact in the limit of large dimensions. In finite dimensions it provides a useful mean-field 
description of this physically relevant phenomenon. 

Two of us (A.S. and P.K.) are grateful to L. Borkowski for stimulating discussions at 
early stages of this work. D.V. acknowledges a useful discussion with P.Horsch. This work 
was supported by the NHMFL (A.S. and P.K.), DOE DE-FG05-91ER45462 (P.K.) and the 
Deutsche Forschungsgmeinschaft through Sonderforschungsbereich 341 (R.S. and D.V.). 

APPENDIX: 

In this appendix we show that our generalization to finite dimensions, Eq. (p4D, retains all 
analytical properties in the upper and lower half planes, as required of Gij{z). Our strategy 
is as follows: First we show that the imaginary parts of ^{z) and z have opposite signs, 
which implies that z and z — S(z) always lie in the same half plane. In particular, z — J]{z) 
has a nonvanishing imaginary part so long as z is not purely real. Then, making use of the 
known analytical structure of Gf^\z) we conclude that: (1) Gij{z) = Gf^\z — ^{z)) has no 
poles in either the upper or lower half planes; (2) '^mG{z) has a fixed sign throughout each of 
the two half planes. The latter point guarantees that the spectral part of G is non-negative. 
Finally, the spectral sum rule (i.e., the normalization of the spectral function to unity) is a 
direct consequence of the fact that G{z) falls like 1/z for large \z\ and has no poles outside 
the real axis. 

We now construct the first and principal step in the procedure outlined above. From the 
asymptotic behavior of Gf^\z) it follows that G{z) = G^'^^{z — Tj{z)) falls like 1/z for large 
\z\, while Gi{z) and G2{z) decay according to higher powers of 1/z. Thus by virtue of Eq. 
(|I^), S(z) ~ ^ for \z\ — s> oo. Asymptotically it is therefore apparent that the imaginary 



18 



parts of ^{z) and z acquire opposite signs. This relation can now be extended to the entire 
upper and lower half planes using the recursion relation, Eq. (|18|). Replacing G{z),Gi{z) 
and 6*2 (-z) with explicit fc-space integrations, Eq. (0) is rewritten as 

1=1 \ ' K ^_ n=f:m n=i 

(Al) 



Here the different terms in the curly brackets correspond to G{z),Gi{z) and G2{z), respec- 
tively. These may be easily brought to the form 
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7* 2f*2 /-TT d. jr. fy;f_ COsf/Crw , 

which clarifies the positivity of the numerator in the integrand. Explicitly taking the imagi- 
nary part of the right-hand side in Eq. ( [A2[ ), we find that QmT.(z + ^) and 53m {z — S(z)} 
always have opposite signs. Thus, given that asymptotically QmJ]{z) is negative (positive) 
in the upper (lower) half plane, and since any point in the half plane can be approached from 
the asymptotic regime by way of repeated applications of the recursion relation, 53mS(2;) 
remains negative (positive) throughout the entire half plane. Hence the imaginary parts of 
z and Ti{z) always acquire opposite signs. 

Having established this point, the analytical properties of Gij{z) follow along the lines 
illustrated above. Both the positivity of the spectral function and the spectral sum rule 
have also been confirmed numerically for several cases, including those displayed in Figs. 7, 
8 and 9. 

Finally, we wish to elucidate the importance of retaining G2{z) in Eq. (|2^). If G2 is 
eliminated from Eq. (p^ we violate, for any finite d, the positivity of the curly brackets in 
Eq. ( [All ). Consequently 53mS(2;) changes signs both within the upper and lower half planes, 
as does the imaginary part of z — S(2;). The branch cut G\^\z) possesses on the real axis 
thus shifts into the upper and lower half planes, and Gij{z) is no longer analytic. 
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FIGURES 

FIG. 1. Hole trajectories that become important in the presence of a ti term, (a) Diagonal 
motion; (b) Triangular loop motion. 

FIG. 2. A schematic description of the paths that contribute to the off-diagonal self-energy, 
Tiij. The hole starts its motion at site j. It approachs the intermediate site / for the first time 
(segment 1) before revisiting site j (segment 2). Then it completes its motion to site i by way of 
/ (segments 3 and 4). 

FIG. 3. Bare DOS due to pure ti-hopping on a hypercubic lattice in d = 2,3, oo, and on a 
Bethe lattice with an infinite coordination number, Z = oo. 

FIG. 4. DOS for J*/t* = 0.4 and tl/t* = 0.05, 0.2 and 0.4 for a hypercubic lattice at d = oo. 

FIG. 5. Similar to Fig. 4, but with tl/t* = —0.1 and —0.2. Notice the disappearance of all 
discrete (resonant) structures for tl/t* = —0.2 . 

FIG. 6. Bethe lattice with coordination number Z = oo : DOS for J*/t* = 0.4 and tl/t* = ±0.4. 
FIG. 7. DOS in d = 3 [see Eq. (H)] : DOS for J* /t* = 0.4 and tl/t* = ±0.2, ±0.4 . 

FIG. 8. Same as figure 7 with J* /t* = 0.1 and tl/t* = -0.4 . 

FIG. 9. The A{k,Lo) spectral function for d = 2 and k = (f,f). Here J* /t* = 0.2 and 
tl/t* = -0.6 . 

FIG. 10. The quasiparticle dispersion, E{k), for d = 2, (a) as a function of J* for fixed 
tl/t* = -0.6, (b) as a function of for fixed J*/t* = 0.1. Here F = (0,0), M = {tt,0) and 
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